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Abstract: We compared two digital elevation models (DEMs) for the Echo Mountain SE quadrangle in the Cascade 
Mountams of Oregon Comparisons were made between 7.5-minute (l:24,000-sca!e) and q 1-degrle (l:250,000-scale) 
presented S VanableS of elevatl0n ' slo P e as P ect ' and ^ gradient. Both visual and statistical differences are 


INTRODUCTION 

T he methods required to assess the extent and distribution 
of natural resources often require that topographic factors 
be taken into account. This may be because these factors are 
directly influencing the resources of interest (vegetation pat- 
terns), or because topographic factors influence the data from 
Earth resources satellites. 

In forested environments, for example, slope gradient and 
slope aspect are known to influence distribution of various plant 
species. The shading resulting from various slope and aspect 
combinations along with solar angle of incidence also influences 
the appearance of forest types on remotely sensed imagery. By 
accounting for topographic variables such as slope gradient, 
slope aspect, and elevation, we should be able to improve our 
capacity to conduct inventories of natural resources. Eby's (1987) 
study of the effect of sun incidence angle data, calculated from 
digital elevation data, is an example in the use of digital Landsat 
data. Stoszek (1977) and Miller and Heller (1978) both deter- 
mined that topographic variables -including elevation, aspect 
and locations near ridgetops- could be used to predict infes- 
tations by the Douglas-fir tussock moth. 

With Geographic Information Systems (GIS), topographic var- 
iables such as these can be developed from Digital Elevation 
Models (DEM) available in digital formats from the United States 
Geological Survey (USGS). Greenlee (1987) developed a method 
to extract drainage networks and other rasterized linework from 
a DEM for input into a vector GIS. Jenson and Dominique (1988) 
reported that DEMs can be useful for hydrologic applications 
such as delineating watersheds and developing overland flow 
path models. Their results indicated that the accuracy and detail 
of the hydrologic information was dependent on the resolution 
and quality of the DEM. 

Two commonly used DEM formats are l:250,000-scale (1-de- 
gree) Defense Mapping Agency (DMA) data and l:24,000-scale 
(7.5-minute) DEM data. The DMA data are available for the whole 
of the conterminous United States, but DEM data are available 
only for selected areas where the USGS has collected data for 
developing and updating their 7.5-minute quadrangle sheets 
(U.S. Geological Survey, 1987). The 7.5-minute DEM data are 
more highly resolved (30 by 30 m per cell) than the 1 -degree 
DMA data, at three arc-seconds (about 65 by 92 m for this mid- 
latitude area), but its limited availability invites the question of 
whether the resolution and accuracy of the DMA data, which is 
available over the entire United States, is adequate for inclusion 
with large area surveys incorporating satellite-acquired data. 

We compared these two data forms on a USGS quadrangle- 
sized area in the Oregon Cascade Mountains. The study area 
chosen was the area covered by the Echo Mountain SE quad- 
rangle (also referred to as Tamolitch Falls) within the McKenzie 
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River drainage east of Eugene, Oregon. The western two-thirds 
of the quadrangle is steep and deeply dissected, and the eastern 
portion is flatter. This area is typical of many forested moun- 
tains of the Pacific Northwest. The objective was to quantify 
the differences between slope gradient, slope aspect, and ele- 
vation images derived from the two digital elevation data forms. 

METHODS 

The 7.5-minute data for the Echo Mountain SE quadrangle 
were read as elevations representing distance above mean sea 
level (MSL) in one-metre intervals with a minimum of 581 m and 
a maximum of 1747 m. The data were displayed as a level-sliced 
color image, with red colors as higher elevations and magenta 
for lower elevations (Plate 1). 

Derived images of slope gradient and slope aspect were com- 
puted on the original 16-bit data (L.N. Brantley, personal comm.; 
see the Appendix for algorithms used). Eight sector aspects and 
a flat category were computed, and slopes were calculated in 
percent. Level-slicing was used to develop a color scheme for 
image representation of these two data sets. 

The 1 -degree data, a full scene of DMA data covering the area 
within 122 and 123 degrees west longitude and 44 and 45 north 
latitude, were read from tape to a disk file. A 151- by 151-pixel 
subset of the data covering the same area as in the 7.5-minute 
data was selected from the full scene with a 601 -m minimum 
and a 1707-m maximum elevation. 

Eighty observations, one from each of the two different ele- 
vation data sets, were obtained from a systematic sampling of 
the study area at points with a 1000-m spacing. A simple linear 
regression analysis was conducted to determine the relationship 
between the 80 elevations from the 7.5-minute data set and the 
1-degree data set. 

For the 1 -degree data, slope gradient and slope aspect images 
were derived, using the same methods described above before 
they were rescaled, rectified, and resampled to 30- by 30-m 
pixels. The resampling was conducted using a nearest neighbor 
routine which used the intensity of the closest pixel to deter- 
mine the value of the output pixel. Slope gradient and slope 
aspect images were colored using the same scheme as with the 
7.5-minute data. 

A difference image was generated by subtracting elevations 
at each pixel location for the whole image area. Comparisons 
for both slope aspect and slope gradient variables were accom- 
plished with an overlay analysis to generate a matrix of obser- 
vations showing frequencies of coincidence of all classes. 

RESULTS AND DISCUSSION 

The relationship between the elevations sampled from the 
7.5-minute and the 1-degree data sets (before resampling) is 
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1 -Degree Data 

■ Warm colors 

are higher 
in elevation 
than cool 
colors. 

Plate 1. Comparison images of 7.5-minute and 1 -degree digital elevation 
Mountain SE Oregon Quadrangle. 


7.5 Minute Data 



models for the variable of elevation for the Echo 




shown in Figure 1. The relationship between these variables 
was found to be linear with a coefficient of determination of 
84.3 percent. The slope, 0.99, and intercept, 36.9, were not sig- 
nificantly different from 1.0 and 0.0, respectively (P < 0.01), in- 
dicating that there was no significant bias. The standard error 
was 30.6 m, and the residuals from the least-squares regression 
line appeared to be distributed randomly. By overlaying the 
elevation layers from the two different data sources, the mean 
elevation difference, on a pixel by pixel basis for the entire area, 
was calculated to be 31.0 m, with a standard deviation of 6.99 
m. Most of the higher differences appeared to be associated 
with steeper slopes. A visual comparison of the 1-degree and 
7.5-minute elevation maps suggest considerable correspon- 
dence of general patterns, but a greater amount of detail in the 
7.5-minute image (Plate 1). 

The visual correspondence of general pattern between the 1- 
degree and the 7.5-minute data sets was also apparent for both 
slope aspect and slope gradient images (Plates 2 and 3). Pattern 
heterogeneity was higher in the 7.5-minute data sets, probably 
due to the finer original pixel resolution. While the modes for 
aspect class comparisons occurred as expected (East with East, 
Southeast with Southeast, etc.) only 36 percent of pixels were 
in the same class in both images (Table 1). Sixtv-seven percent 
of pixels were in the same or adjacent (East with Northeast, 
East and Southeast, e.g.) classes. In the slope gradient images, 
29 percent of all pixels were in the same class in both images 
(Table 2), and gradients calculated from the 1-degree data were 
lower than those from the 7.5-minute data. The differences ap- 
parent in the images seemed to derive from a greater level of 
detail in the 7.5-minute than in the 1-degree data. 

While we have noted anomalies in 1-degree data in other 
areas, of low relief, which create obvious artifacts in derived 
topographic images, no instance of this w r as found in the data 
used in this comparison in areas of mountainous terrain. This 
is true also of other 1-degree data sets we have used. For ap- 
plications in mountainous terrain, 1-degree data seem an ap- 
propriate match for Landsat Multispectral Scanner data and for 
databases where data storage limitations may be of concern. 



Elevations from 7.5 minute data 
(meters X 100) 

Fig. 1 . Scatter diagram of elevations obtained from 
a 7.5-minute digital elevation model (x axis) and a 
1 -degree digital terrain model (y axis) on the Echo 
Mountain SE quadrangle area of the central Ore- 
gon Cascade mountains. 


One final note has to do with the order of processing for 
digital terrain data sets. Resampling to adjust pixel size, recti- 
fication of images, and rescaling to different intervals resulted 
in different final image values depending on the order of 
processing. In some cases the differences were great. In this set 
of analyses, for example, we retained the 16-bit format as far 
along in the chain of processing as possible, and final images 
for slope and aspect yielded very different quantitative esti- 
mates of flat terrain. Because the original data are in one-metre 
intervals, and the algorithm for slope aspect determines a non- 
flat aspect for any unequal elevation, very few slope aspect flat 
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Table 1. A Comparison of Slope Aspects Derived from 7.5-Minute and 1 -Degree Digital Elevation Models (DEM). Numbers Are 
Counts of Coincident Occurrence for the Entire Echo Mountain SE Quadrangle. 
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34% 

46% 
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Table 2. A Comparison of Slope Gradient Classes Derived from 7.5-Minute and 1 -Degree Digital Elevation Models (DEM). Numbers 
Are Counts of Coincident Occurrence over the Entire Echo Mountain SE Quadrangle 
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pixels resulted. In slope gradient determinations, rounding can 
result in zero percent calls without absolute equality in the el- 
evations used in calculations, and the zero percent classes were 
greater: 10,715 zero percent slopes versus 4161 "flat" aspects in 
the 1-degree data, and 4,799 zero percent slopes versus 425 "flat" 
aspects in the 7.5-minute data. Users of digital terrain data should 
carefully consider and check each step in their processing to 
ensure meaningful results. The methods of production can ef- 
fect the kinds of artifacts that occur and the accuracy of the 
results. 

We have concluded that, in our applications in wildlife hab- 
itat analysis, 1 -degree DEM data are preferable to our developing 
our own digital elevation data for the extensive areas for which 
7.5-minute data are not available. We suggest, however, that 
additional comparative research be conducted involving both 
types of DEMs for a variety of terrain types. The results reported 
in this paper may or may not be indicative of DEM differences 
found in other areas, but the comparison techniques in this 
paper can be applied to other areas with different topographic 
profiles. 
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APPENDIX 

Computations of slope gradient and slopect aspect (source, 
L.N. Brantley, personal comm.): 

consider the 3- by 3-pixel matrix of elevations around the point 

fry)- 
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Computation of average x and y elevation changes; 


Ax l = c- a 
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^ _ Ax l + Ax 2 + Ax 3 
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(1) Computation of slope gradient: 
let s = cell size 


Ae = \/Ax 2 T Ay 2 

if Ae<s, percent slope = Ae/s* 100 

if Ae>s, percent slope = 100 + [(1 - s/Ae)* 100] 

(2) Computation of slope aspect: 

if Ax = 0 and Ay = 0, aspect = FLAT, 
otherwise 


6 = tan - 
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FOUNDATION 


Dedicated to the improvement of the quality of life 
on earth through better use of geographic information. 


The International Geographic Information Foundation was originally founded in 1979 by the American Society of 
Photogrammetry to administer the funding of educational scholarship and awards programs. In 1986, its scope was expanded 
in response to the need by other professional societies concerned with geographic information for assistance and expertise in 
r aiding and administering funding for scholarship, awards, and research programs. 


The IGI Foundation exists to: 


• Conduct programs in the field of the mapping 
sciences which provide an avenue for practitioners 
of the geographic sciences to make tangible contri- 
butions from their own knowledge and experience to 
the advancement and practices of the sciences. 

• Extend to other organizations the opportunity to 
exchange information and mapping experience, and 
the practical solutions of mapping sciences 
problems. 


• Encourage improvements in the geographic 
information sciences and their applications by 
making grants for the conduct of research in the 
fields of the mapping sciences and establis hing 
awards programs to recognize significant 
contributions to humanity through the sciences 

• Improve public awareness of the benefits available 
from the geographic information sciences and of the 
contributions they make to society. 


The Foundation aids the societies by assisting in raising funds through charitable contributions and endowments for 
worthwhile educational, research, and awards programs, and providing expertise in funds management. The Foundation solicits 
charitable contributions from individuals, corporations, interested organizations, and other foundations. Each society served 
by the Foundation develops its own programs and performs its own internal selection of scholarship, awards, and research gran 

recipients. 


A variety of contribution options are available: 


• Established Fund Awards - Named Awards 
Programs may be established and maintained 
perpetually with a minimum donation of $15,000, 
which may be contributed over a 3-year period, or 
the fund will be disbursed as a one-time award. 

• Birdseye Club - Named in honor of Colonel 
Claude B. Birdseye, a distinguished mapping 
scientist and topographer, and established to honor 
major individual and corporate benefactors of the 
Foundation (Individual contribution-$ 1,000; 
Corporate contribution-$5,000). 

• General Contributions - Accepted with 
appreciation. 


• Centurion Club - Established to recognize 
individual benefactors who make an annual, 
unrestricted contribution to the Foundation and 
corporate benefactors who donate a premium to the 
Foundation’s benefit auction (Annual unrestricted 
Individual donation of $25, $50, or $100; Annual 
unrestricted Corporate donation of premium to 
Benefit Auction having a minimum retail value of 
$250 or 500). 

• Gift Annuity Program - Charitable 
Remainder Annuity Trust for the benefit of the 
foundation (a deferred giving program). 


To learn more about the Foundation, how to establish an awards, scholarship, or research program and how you can 
contribute, contact: 


• International Geographic Information Foundation • 
5410 Grosvenor Lane, #210 
Bethesda, MD 20814-2160 
301-493-0290 
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